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Abstract Hyper-velocity stars are believed to be ejected out from the Galactic center through 
dynamical interactions between (binary) stars and the central massive black hole(s). In this 
paper, we report 19 low mass F/G/K type hyper-velocity star candidates from over one mil¬ 
lion stars of the first data release of the LAMOST general survey. We determine the unbound 
probability for each candidate using a Monte-Carlo simulation by assuming a non-Gaussian 
proper-motion error distribution, Gaussian heliocentric distance and radial velocity error dis¬ 
tributions. The simulation results show that all the candidates have unbound possibilities 
over 50% as expected, and one of them may even exceed escape velocity with over 90% 
probability. In addition, we compare the metallicities of our candidates with the metallicity 
distribution functions of the Galactic bulge, disk, halo and globular cluster, and conclude that 
the Galactic bulge or disk is likely the birth place for our candidates. 


Key words: stars; low-mass — stars; kinematics and dynamics — Galaxy; abundances — 
stars; fundamental parameters — stars; distances 


1 INTRODUCTION 


The hypervelocit y stars 


(HVS s), discoyered in the Galactic halo, are trayelling so fast that they can escape 


from the Galaxy. iHillsI (1198811 first predicted their existence, and the existence are the powerful eyidence 


of massiye black hole at the center of the Milky Way. A natural explanation is that they may be ejected 
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out from the Galactic center (GC) by int eractions of stars with the massive black hole (MBH) or the hy¬ 
pothetical binary MBHs as predicted by 


Hillsl (119881) and 


Yu & Tremainel (120031) . Such ejection mecha¬ 


nism^ canbedividedJnto_thre£_categories^_tidalbre^up of binary stars in the vicinity of a single MBH 


dHills 


1988 


Yu & Tremaine 


2003 


Bromley et al. 


20061) . and the bi nary stars are 


the vicinity of the MBH from the young stellar disk in the GC (e.g., 


from the Galactic bulge (IPereti 


Sesana et al.ll2007 


Loo?: Merritl 




Lu et al 


probably injected i nto 


201C; 


[g); single star encounters with a binary MBH (lYu & Tremainell2003 


Zhang et al 


2010 or 


inJ? 


MemtgHOOw or single star encounters with a cluster of stellar mass black holes around 


the MBH (lO’Learv & Loebll2008h . 
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However, the black hole acceleration mechanism cannot explain a type of HVSs such as US 708 (Hirsch 
et al. 2005), which is not originated from the center of our Galaxy (e.g., Geier et al. 2015). This type of 
HVSs is likely to be the ejected donor remnant of a thermonuclear supernova in the white dwarf + helium 
star scenario (see Wang & Han 2009; Geier et al. 2015), in which LAMOST are searching for such type of 
HVSs based on the theoretical results of Wang & Han (2009). 

Alternatively, other ejection models can also accelerate st ars to h i gh speed. For example, binary 


disruption in the dense s tellar clusters such as Galactic disk (IBlaauw 


1961 


Leonard & Dewev 


Naniwotzki & Silva 


1993 


2012h . in such case, a supernova explosion of more massive evolved component can 


accelerate its companion to high speed. Tidal disruptions of dwarf galaxies in th e Milky Way can also 


produce high velocity stars (lAbadi et al. 


2009; 


Lu et al 


201C; 


Tevssier et al 


20091) . such mechanism can 


produce high speed star stream or a population of old isolated ‘escaped’ (unbound) or ‘wondering’ (bound) 
stars. 


Seventeen years after Hill’s predictio n, three HVSs were successively discovered (IBrown et al. 


2005 


Hirsch et al 


2005 


Edelmann et al 


20051) . they are m assive O or B type stars located in the Gala ctic halo. 


Until recently, over 20 unbound HVS were identified (IBrown et al.ll2009 


2012 : 


Zheng et al.ll2014l) . most of 


them are massive B type stars, and are located at distant Galactic halo with Galactocentric distances larger 


than 25 kpc. An interesting exception is a identified B type HVS discovered bv IZheng et al 


(I2OI4I) . which 


is the first HVS discovered in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) 
survey, and it is the brightest HVS currently known, and locate at a Galatocentric distance of 13 kpc. 

Assuming a Salpeter initial mass f unction (IMF), the expected solar mass H VSs are about 10 times 


more abundant than the 3^ Mq HVSs dBrown et al 


20091) . 


Kollmeier et al. 


(l2009t) systematically searched 


for such low mass HVSs in about 290,000 spectra of the Sloan Digital Sky Survey (SDSS) (York et al. 
2000; Gunn et al. 1 998, 2006), ho wever, they found only 6 metal-poor stars that can be possibly taken 


as HVS candidates. 


Li et al. 


(I2OI2I) also searched F and G type HVSs from over 370,000 stellar spectra 


of data release seven of the SDSS, and they presented a low mass metal-poor HVS candidate catalogue. 


Palladino et al 


(I2OI4I) identified 20 G and K type HVS candidates from approximately 240,000 stars of the 


Sloan Extension f or Galactic Underst anding and Exploration (SEGUE) (Yanny et al. 2009) G and K dwarf 
samples. Besides, 


Zhong et al. 


(I2014h reported a catalog of 28 high velocity star candidates from the first 


data release (DRl) of the LAMOST general survey (Luo et al. 2015), which cover a much broader color 
range than ever, and 17 of them are F, G, or K type low mass stars. These current searching results might 
suggest that the IMF of the parent population of these HVSs is top heavy. A top heavy IMF of the HVS 
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parent population is 


Kollmeier et al 


possibly consistent with the disk origination (IBartko et al.ll2010t IZhang et alJ 12013 


20101) . In order to distinguish the ejection mechanisms of HVSs and put constraints on the 


origin of the parent population of HVSs, it is quite necessary to search for the low-mass HVSs. 

In this paper, we systematically search and investigate HVSs with stellar spectra of LAMOST DRl, 
and totally find 19 HVS candidates. The structure of this paper is as follows. In Section 2, we introduce the 
LAMOST and DRl data productions in detail. In Section 3, we present a series of spectroscopic, photomet¬ 
ric and dynamic criteria to select HVS candidates. In Section 4, we analyze the probability for each HVS 
candidate that they can escape from our Galaxy. In Section 5, we compare the metallicities of our 19 HVS 
candidates with the [Fe/H] distributions of Galactic bulge, disk, halo and globular cluster, and conclude 
that our HVS candidates are likely originated from the Galactic bulge or disk. Finally, the discussion and 
conclusion are given in Section 6. 

2 THE FIRST DATA RELEASE OF THE LAMOST GENERAL SURVEY 


The LAMOST is a 4 meter quasi-meridian reflecting Schmidt telescope, it adopts novel active optics tech¬ 
nique, which allows both a large effective aperture of about 4 m and a wide held of view of 5°. The focal 
surface of LAMOST has 4000 precisely positioned optical fibers, which are equally connected to 16 spec¬ 
trographs, thus it can observe 4000 targets simultaneously. Each spectrograph is equipped with a red channel 
CCD camer a and a blue one, which can simu l taneously provid e red and blue spectra of observed targets 


respectively (ICui et al.ll2012 ; 


Wang et al.lll996t 


Su & Cuill2004 . 


The primary scientific goal of LAMOST survey is to investigate the large-scale structure of the universe, 
as well as structure and evolution history of Galaxy, and it consists of two main parts. The first part is the 
LAMOST Extra-Galactic Survey (LEGAS) of galaxies, and the second part is t he LAMOST Experiment for 


Galac tic Understanding and Exploration (LEGUE) Survey of the Milky Way (IZhao et al.ll2012 ; 


Deng et al 


20121) . Considering the science goals and the targets available, the LEGAS is consists of galaxy survey and 


QSO survey, and the LEGU E is divided into three parts, i.e., t he Galactic antice nter survey, the disk sur¬ 


vey, and the spheroid survey (ILiu et al 


2014 


Chen et al.ll2012: 


Carlin et al. 


2012). The Galactic anticenter 


survey (LSS-GAC) covers a significant volume of the Galactic thin/thick disks and halo in a continuous 
sky area of ^ 3,400 square degree, which is centered on Galactic anti-cente r area with Galactic longitudes 


150° < 1 < 210° and latitudes |6| < 30° ( Liu et al. 


2014 


Yuan et al 


20141) . The disk survey choose eight 


low and bright plates along t he Galactic plane which are nearly uniformed distributed in the region 0° < 

len et alJl^ 


which plan to obs erve 5.8 million objects with r < 16.8 at |5| > 30° (lYang et al. 


Zhang et al.ll20O . 


Yang et al. 

2012 

Deng et al. 

2012 


After one year pilot survey (ILuo et al 


2 OI 2 I) . LAMOST began its first year general survey from 


September 2012, and completed its first data release (DRl) to domestic data users and foreign collabo¬ 
rators in June 2013. DRl totally published 2,204,696 wavelength-calibrated and relative flux-calibrated 
spectra, which are consists of 1,944,329 stars, 12,082 galaxies, 5,017 quasars, and 243,268 unknown ob¬ 
jects. These spectra cover a wavelength range of 3690—9100 A with a resolution of R ^ 1800. In addition, 
DRl also published five spectroscopic parameter catalogs, which are the General catalog, the A, E, G and K 
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Fig. 1 The Hertzsprung-Russell diagram of F, G and K type stars of LAMOST DRl, where blue, 
green and red points represent F, G, and K type stars respectively. The black dot-dash line shows 
the location of log^ = 3.75 for F type stars, and the black dashed line displays logf; = 4.0 for G 
and K type stars. 


type stars catalog, the A type stars catalog, the M dwarfs catalog and the observed plate information catalog 
respectively. The A, F, G and K type stars catalog provides effective temperatures (T eff), surface gravi¬ 
ties (logp), metallicities ([Fe/H]) and heliocentric radial velocities for 1,061, 918 stars, these s pectroscopic 


parameters are indispensable for searching and studying hyper-velocity stars (ILuo et al.ll2015l) . 


3 IDENTIFICATION OF THE LAMOST DRI HVS CANDIDATES 

Our HVS candidates are drawn from F, G and K dwarfs of the LAMOST DRl, we use five steps to search 
for them, and the numbers in below brackets indicate the number of stars left after each step. 


1. Selecting F, G and K dwarfs from LAMOST DR1[519,027]. The LAMOST DRl officially released 
atmospheric parameters for 1,061, 918 A, F, G and K type stars, which are derived by the LAMOST 


Stellar Parameter Pipeline (LASP) (ILuo et al.ll20L5h . We hst select F, G and K type stars with 3600 < 


Teff < 7500, and plot their Hertzsprung-Russell diagram in Figure [T] From Figure [T] we can see that 
F, G and K dwarfs can be initially selected with the criteria: 1) logp > 4.0, when 3600 < Teff < 6000; 
2) logp > 3.75, when 6000 < Teff < 7500, and we totally obtain 737,023 F, G and K type dwarfs. 
Then, we upload equatorial coordinates of our 737,023 dwarfs to the ‘MyDB’ database of SDSS DRIO, 
and obtain photometry information of 521,618 dwarfs from the ‘PhotoObjalT table, as well as proper 
motions of 519,027 dwarfs from the ‘ProperMotions’ table. Finally, 519,027 F, G and K dwarfs with 
both photometries and proper motions are selected. 









Hyper-velocity star candidates from LAMOST 


5 


2. Further selecting F, G and K dwarfs with the color and magnitude methods[ 191,405]. 

(a) go < 20.2, tq < 19.7, 0.2 < {g — r)o < 0.48, forF dwarfs 
ro < 19.7, 0.48 < (g — r)o < 0.55, for G dwarfs 

ro < 19.0, 0.55 < (g — r)o < 0.75, for K dwarfs 

(b) A(r) < 0.5 mag 

(c) \b\ > 10 

(d) psfMagErr_g/r/i < 0.05 mag, mode = 1 and dean = 1 

(e) 0.2 mag < {g — i)o < 4.0 mag 


Newby et al 


j201 ih investigated the photometric uncertainties of SDSS, and pointed out that they are 


constant u p to r-band apparent mag nitude (rg) of 19.7, thus we use above criterion (a) to select F, G and 


K dwarfs. 


Schlesinger et al 


(1201211 mentioned that above and below the Galactic plane undergo small 


amounts of extinction, this small amount of reddening can affect target selection, so we use criteria (b) 
to retain F, G, and K dwarfs samples with extinction in r-band (A(r)) less than 0.5 mag (A(r) < 0.5 
mag). As we known that Schlegel et al. (1998) map has a limited spatial resolution and fails at low 
Galactic latitudes, extinction from the Schlegel et al. (1998) map may not represent the true value of 
targets which may lead to inaccurate distance estimates, thus we use criterion (c) to remove dwarfs with 
|6| <10. In addition, in order to make sure SDSS photometry reliable, we use criterion (d) to constrain 
SDSS psf mag errors in g, r, and i band lower than 0.05 mag, and the value of photometric flag ‘mode’ 
and ‘clean’ are 1. In the end, we constrain our dwarf samples to the color range of 0. 2 < (g-'i)o <4.0 
m ag, which is the rel iable range to estimate absolute magnitude with the photometric parallax relation 


of 


Ivezic et al 


(l2008ll . 


3. Calculate phase space coordinates and escape velocities[191,405]. We first calculate heliocentric dis¬ 
tances for 191,405 F, G, and K dwarfs using the distance modulus; 

2Q0.2x(ro-Mr) 


c?(kpc) = 


100 


( 1 ) 


where d is the heliocentric distance in unit of kpc, rg is the r-band deredden apparent magnitude directly 
from the SDSS DRIO ‘PhotoObjAll’ table, and M^. is the r-band absolute magnitude. 

The above step constrains our dwarfs in the color range 0.2 < {g — i)g < 4.0 mag, thus we use the 


photometric parallax relation proposed by 


Ivezic et al 


(120081) to estimate M^. 


Mr{{g - i)o, [Fe/H]) = M°{{g - i)o) + AMr{[Fe/H]) 


( 2 ) 


where M^{g — i) is the color-magnitude relation, AMr([Fe/Hl ) is the absolute ma gnitude correction. 


and they can be calculated by the equations (A7) and (A2) of 


Ivezic et al 


(l2008h . respectively. The 


(ff “ Oo color can be available from the table ’PhotoObjAll’ of SDSS DRIO, and the metal abundance 
[Fe/iT] can be obtained from the A, F, G and K type stars catalog of LAMOST DRl. 

Then, we use a right-handed cartesian coordinate system centered on the GC to calculate Galactic three 
dimension (3D) space positions, the X axis pointing from the Sun to the Galactic center (GC) with 
the Sun at a; = -8 kpc, the Y axis pointing in the direction of rotation and the Z axis pointing towards 
the Northern Galactic Pole. We use a similar coordinate system to calculate Galactic 3D velocity, and 
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assume that the motion of the local standard of rest (LSR) is 220 km s~^, and the velocity of the Sun 
with respect to the LSR is (11.1 km s“^, 12.24 km s“^, 7.25 km s“^) ( Schonrich et akEoiO ). 

To determine which dwarfs are unbound to the Galaxy, we adopt five different Galacti c potential mod 


hereafter XueOS), a spherically symmetric divergent model ( 

Kenvon et a 

. 2008. hereafter KenvonOS), 

two axisvmmetric divergent models (Paczvnski 199C; Konosovetal. 

201ol hereafter Paczvnski90 

and KoposovlO), and a triaxial convergent model (Gnedin etal. 2005. 

hereafter Gnedin05). Beside 

nent bulge-disk-halo models. 

a simple model of KenyonOS, other four potentials are all three-compo 


XueOS, KoposovlO and Gn edin05 adopt a spherical b ulge (IHernauist 


the Miyamoto-Nagai bulge (IMivamoto & jjag ai 


GnedinOS use the Miyamoto-Nagai disk! 


19901), which is different from 


19751) of Paczvns ki90. Paczynski90, KoposovlO, and 


Miyamoto & Nagailll975l) . while XueOS adopt an exponential 


disk. For the halo component, the four thre e component models are entirely different (iNavarro et al. 

][ 


1996t 


Paczvnskill 19901: 


Fellhauer et alJl2006h . For the convergent Galactic potential model, there exists 


i true escape velocity, and we use equation[3to estimate 14sc; 


|t^esc(r)| = ^2 X |4>(r)| 


(3) 


However, for the divergent potential, any HVS with finite space velocity could not really escape from 
the Galaxy. We thus def ine unbound stars as ones which can reach r > 200 kpc with v > +200 km s“^ 
as shown in eauationl4l( Kenvon et al.ll^O^ ): 


|ue.c(r)| = xl2 X (- X 2002 + ||$(r)| - |<h(200)||) 


(4) 


4. Find HVS candidates with clean or reliable proper motions [32]. In step two, we use a series of criteria to 
make sure that our dwarf samples have reliable SDSS photometric parameters, which affect heliocentric 
distance estimate accuracy. In addition, the proper motion distribution of SDSS tends to display large 
proper motion errors, we thus need to ensure that proper motions of dwarf samples are real rather than 
the product of large error s. 


Munn et al 


(12004 . 


2008h presented an improved prop er-motion catalog, which matched each SDSS 


point source to the USNO-B catalog (IMonet et al 


20031) . The SDSS+USNO-B catalog is 90% complete 


to p < 19.7, and has statistical errors of roughly 3-3.5 mas yr ^ an d approximately 0.1 mas yr ^ 


systematic errors for each component proper-motion. 


Munn et al. 


(120041) also defined a series of criteria 


to make sure t hat th e proper-motions from SDSS catal og are reliable, these c riteria were later revised 


by 


Kilic et al. 


(l2006h for their white dwarf samples, and 


Palladino et al. 


(120141) used the revised criteria 


when they search for G and K type HVS candidates from the SEGUE. In this paper, we use the criteria 
to select HVS candidates with ‘clean’ or ‘reliable’ proper-motions, the ‘clean’ proper-motion is defined 
as follows: 

- match=l, which represents only a USNO-B object is within 1” radius of the SDSS target. 

- sigRA < 52 5 and sigDEC < 525, which means that the proper-motion fit must have rms 
residuals less than 525 mas in both right ascension and declination directions. 

- nFit > 6, which shows at least six observations (SDSSh-USNO-B) have been used to determine 


the proper motion. 
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- dist22 > 7, which suggests the distance to the nearest neighbor with g <22 must exceed 7”. 
and the ‘reliable’ proper-motion is dehned as follows: 

- match=l, sigRA < 525 and sigDEC < 525. 

- nFit = 6anddist22 < 7;or, nFit = 5anddist22 > 7 


Kilic et al 


(l2006h pointed out the contamination rate for a target with a ‘reliable’ proper motion is not 


larger than 1.5%. Finally, we totally hnd 32 HVS candidates which can escape from the Galaxy in at 
least one Galactic potential model mentioned in step three, 15 of them have ‘clean’ proper motions, and 
other 17 candidates have ‘reliable’ proper motions. 

5. Finally, HVS candidates are selected with high quality spectra and reliable atmospheric parameters[19]. 
After above four steps, 32 HVS candidates are initially selected from over 190,000 F, G and K dwarf 
samples, we further visually inspect their LAMOST spectra, and hnd ten proper-motion ‘clean’ can¬ 
didates and nine proper-motion ‘reliable’ candidates, which hav e high quality spectra, and their r- 


band signal-to-noise ratio ‘SNR_r’ are listed in Table [T] 


Lee et al. 


(l2008h presents the SEGUE Stellar 


Parameter Pipeline (SSPP), which was used to derive the fundamental stellar atmospheric parameters 
(Teff, logg, and [Ee/H]) for A,E,G and K type stars using multiple methods. We use the version of the 
SSPP used for the seventh data release of the SDSS to verify atmospheric parameters for 19 HVS candi¬ 
dates, and conclude that the parameters derived by the SSPP are roughly consistent with those obtained 
by the LASP. 


Here, we hnally select 19 HVS candidates, all of them have reliable photometry and stellar atmospheric 
parameters, over 98.5% probability of robust proper-motions, high quality spectra and no visual blending. 
Their fundamental parameters such as equatorial coordinates, r-band dereddened apparent magnitudes and 
atmospheric parameters, are shown in Table [T] and their heliocentric distances. Galactic distances. Galactic 
total velocities, and escape velocities obtained by hve Galactic potential models are listed in Tabled 

Prom Table[T] we can see that the value of [Pe/H] errors obtained from the LAMOST parameter catalog 
seem to be large, even much larger than their true value, such large errors will effect the distance error 
estimates. Prom Tabled we can see that the Galactic distance errors are really large, and a fraction of them 
can reach up to or over 10%. Luo et al. (2015) pointed out that external errors released in the LAMOST 
parameters catalog are larger than the real measurement errors, because they rescaled the external errors 
using a ratio. They compared the LASP parameters with high resolution spectra results, SDSS DR9 results, 
and results from Gao et al. (2015) and LSP3(Xiang et al. 2015), the mean external error is 0.125 dex which 
is much close to the true error. Considering this mean error as [Pe/H] errors of our 19 candidates, the 
corresponding Galactic distance errors listed in the brackets of the third column in Tableware really much 
smaller. 

Theoretically, Ivezic et al.(2008) presented out that a 1.0 dex [Pe/H] error will bring about a 1 mag 
absolute magnitude (Mr) error at the median thin-disk metallicity ([Pe/H] = -0.2), and a 0.57 mag Mr error 
at the median halo metallicity ([Pe/H] = -1.50), which represents the mean LAMOST external [Pe/H] error 
of 0.125 dex will result in a Mr error of 0.125 mag at most, and an approximately 6% Galactic distance 


error. 
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Besides, we check whether our candidates are new fin dings.iKollmeier et alJ (120091) present 6 F/G type 


metal-poor HVS candidates from over 290,000 SPS S stars. 


poor unbound HVS candidates from SPSS PR7, and 


Li et al 


Palladino et al 


(1201^) proposed 13 F and G type metal- 


(12014h found 20 G and K type unbound 


HVS candidates in SEGUE G and K dwarf samples from the SPSS PR9. Our 19 HVS c andidates are not in 
the three HVS catalogs by checking equatorial coordinates. Besides, [ 


Zhong et al 


(12014ft presented a catalog 

of 28 high-velocity star candidates from LAMOST PRl, they use a velocity criterion |ru| > 200 km s“^ 
and |l^t| > 300 km s“^ when selecting HVS candidates, where rv is the heliocentric radial velocity, 
and Vgt is the 3P total velocity, and this criterion prevent our candidates exc ept one to be found . Because 


the exceptional candidate is not in the u — g and r — i color ranges where 


Zhong et al 


(120141) estimate 


accurate photometric metallicity, it is also included in our HVS candidates catalog. Among their 28 HVS 
samples, there are 12 stars with spectra type earlier than E type, and other 16 candidates are eliminated by 
our photometry and proper-motion criteria. In summary, we can conclude that our 19 low mass E/G/K type 
HVS candidates are new findings. 


4 ESTIMATING THE RELIABILITY OF OUR CANDIDATES WITH MONTE-CARLO 
METHOD 

Although we use a series of criteria to ensure our HVS candidates with reliable photometry, atmospheric 
parameters and proper-motions in section [2 it is still premature to say that the final HVS samples do 
not contain false-positive detections. With this in mind, we thus consider the probability that our HVS 
candidates are in fact unbound to the Milky Way. To obtain such an unbound probability for each HVS 
sample, we built a Monte-Carlo simulation to sample a million realizations of orbit parameters. 


Dong et al 


(1201 ift present a non-Gaussian probability distribution function (PDE) of proper-motion 


errors using quasar samples with ‘clean’ proper-motions, which contains a Gaussian core and an ex¬ 
tended wing. Applying this error distribution model, we randomly produce a million total proper-motion 
errors (pmerror) with the inverse function method. Assuming the proper-motion error is isotropic in 
the ra-dec plane, we can produce a million angles ‘0’ using an uniform distribution model, and obtain 
a million proper-motion errors in the ra and dec direction using pmraeinor = pwerror x cos{6) and 
pmdecerroi- = pWerror X sin(0), where pmraerroi- and pmdecerroi- are random proper-motion errors in 
ra and dec directions respectively. Using the proper-motion measurements from SDSS in ra and dec direc¬ 
tions, two component proper-motion errors from SDSS and the random proper-motion errors pmraerroi- and 
pmdecerror, we can obtain a million random two components proper-motions. 

In addition, we randomly generate a million radial velocities and heliocentric distances assuming a 
Gaussian error distribution function, and a million random Galactic 6D phase space coordinates and escape 
velocities at each Galactic distance in the million realization can be further obtained. In such a Monte-Carlo 
simulation, unbound probability for each HVS candidate can be derived by the fraction ^o^o(^ooo° ’ where 
Vgt and 14sc are Galactic total space velocities and escape velocities respectively, and A^Vgt>Yeso ’^he 
number of Vgt larger than 14sc in a million realizations. 

Table[3]lists unbound possibilities for each HVS candidate, and ‘—’ in this table represents that the can¬ 
didate is not unbound in this potential model, and we do not calculate unbound possibilities in these cases. 
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Table 1 Fundamental parameters of 19 HVS candidates. 


HVS 

Designation 

ra“ 

(degree) 

dec“ 

(degree) 

rO'’ 

(mag) 

SNRj-= 

tvs'* 
(km s“^) 

flaCOs{6Y 
(mas yr“^) 

^5^ 

(mas yr“^) 

Teff-^ 

(K) 

loggf 

[Fe/H]^ 

[Mg/Fe]9 

1 

J172524.12+565709.6 

261.3505 

56.95267 

13.95 

101 

-103±11 

-17.19±2.62 

88.11±2.62 

5251 ±102 

4.65±0.43 

-0.76±0.36 

0.43 

2 

J170333.23+373102.3 

255.8885 

37.51733 

17.50 

29 

-60±10 

-22.60±3.00 

18.79±3.00 

5069±156 

4.40±0.41 

0.08±0.38 

-0.05 

3 

J132422.30+312841.6 

201.0929 

31.47823 

15.76 

18 

-23±16 

-24.35±2.36 

-27.53±2.36 

5828±330 

4.03±0.66 

0.02±0.62 

0.33 

4 

J091255.48+140413.8 

138.2312 

14.07052 

15.26 

13 

46±22 

9.82±2.52 

-54.49±2.52 

6231 ±307 

4.67±0.36 

-0.86±0.82 

0.57 

5 

J130548.65+282410.7 

196.4527 

28.40298 

17.43 

16 

114±12 

-21.49±2.81 

-37.11±2.81 

5931 ±447 

4.28±0.93 

-1.61±1.536 

1.11 

6 

J133115.50+150438.9 

202.8146 

15.07748 

18.49 

62 

-29±6 

-20.87±3.05 

4.08±3.05 

4854±90 

4.68±0.29 

0.12±0.24 

-0.03 

7 

J175513.55+511927.4 

268.8065 

51.3243 

13.87 

38 

-72±11 

8.06±2.73 

48.98±2.73 

5228±170 

4.45±0.47 

-0.34±0.46 

0.32 

8 

J113116.03+571131.1 

172.8168 

57.19199 

16.19 

23 

-91±11 

-43.52±2.64 

-38.49±2.64 

6083±259 

4.07±0.50 

-1.47±0.89 

0.58 

9 

J121811.06+284659.9 

184.5461 

28.78333 

14.08 

54 

-18±8 

-48.94±2.57 

-0.32±2.57 

5068±117 

4.64±0.35 

0.04±0.30 

-0.03 

10 

J115209.12+120258.0 

178.038 

12.04946 

15.86 

22 

206±15 

-32.11±2.51 

19.24±2.51 

5669±280 

4.11 ±0.64 

-0.01±0.57 

0.10 

11 

J004028.68+393853.0 

10.11953 

39.64808 

16.13 

19 

-54±23 

-34.08±2.57 

-23.19±2.57 

6119±290 

4.31 ±0.50 

-0.51 ±0.69 

0.42 

12 

J171952.43+525035.6 

259.9685 

52.84325 

15.62 

24 

-67±13 

-32.84±3.2 

83.20±3.2 

5703±244 

4.20±0.59 

-0.11 ±0.51 

0.27 

13 

J063934.38+280912.8 

99.89327 

28.15358 

17.02 

8 

15±30 

0.11±2.52 

15.04±2.52 

6064±431 

3.94±0.59 

-0.33±0.88 

-0.18 

14 

J005233.53+413322.6 

13.13972 

41.5563 

14.81 

23 

-30±23 

-48.83±2.55 

-32.97±2.55 

6187±282 

4.34±0.48 

-0.43±0.64 

0.02 

15 

J012947.93-021343.2 

22.44971 

-2.228684 

14.85 

50 

13±13 

-24.22±3.82 

30.44±3.82 

5860±183 

4.25±0.56 

-0.36±0.42 

0.26 

16 

J075303.30+272657.0 

118.2638 

27.44918 

17.81 

14 

75±32 

-18.13±3.00 

4.51±3.00 

6017±399 

4.26±0.70 

-0.39±0.75 

0.36 

17 

J142235.20+455631.3 

215.6467 

45.94204 

14.89 

44 

-121±12 

-42.77±2.43 

11.14±2.43 

5231±138 

4.68±0.43 

-0.57±0.44 

0.38 

18 

J130744.34-004449.5 

196.9348 

-0.747102 

17.18 

9 

14±20 

-9.72±3.11 

11.48±3.11 

6035±408 

4.28±0.50 

-0.27±0.86 

0.56 

19 

J103858.44+565558.1 

159.7435 

56.93283 

13.74 

74 

-15±14 

-88.10±5.61 

97.36±5.61 

5785±146 

4.28±0.54 

-0.74±0.41 

0.36 


Notes. The candidates HVSl—HVSIO all have ‘clean’ proper-motions, and the candidates HVSll— 
HVS19have ‘reliable’ proper-motions. 

“ Equatorial coordinate from the SDSS ‘PhotoObjAlT catalog. 

Dereddened r band apparent magnitude from the SDSS ‘PhotoObjAlT catalog, 
r-band signal to noise ratio from LAMOST parameter catalog. 

Heliocentric radial velocity from LAMOST pai'ameter catalog. 

® Proper motion in both right ascension and declination directions from the SDSS ‘ProperMotions’ catalog 
■f Atmospheric parameters from the LAMOST parameter catalog. 

3 Note that there are large uncertainties in [Mg/Fe] measurements. 


From this table, we can see that unbound probability for each HVS candidate exceeds 0.5 as expected, and 
it is even over 0.9 for HVSl in the XueOS potential model. In addition, for each HVS candidate, the value 
of unbound probability depends on the adopted Galaxy potential model to a certain extent, such as HVSl, 
the probability varies from 0.64 in the GnedinOS model to 0.93 in the XueOS model, and the probability 
difference between two potential models changes from 0.07 to 0.29. Among our 19 HVS samples, only 
seven candidates are unbound in all five potential models. Figure |2] shows their distribution of Vga,i — 14sc, 
where V^ai and Ksc are total Galactic velocities and escape velocities of each realization. From this figure, 
we can see that the total velocity exceeds escape velocity in most cases. 

Actually, the value of unbound probability mainly relies on the difference between total space velocity 
and escape velocity. When the total velocity is much larger than the escape velocity, the effect of parameter 
error and Galactic potential model will be extremely small. Conversely, when total velocity is just larger 
than the escape velocity, parameter error and potential model will greatly affect the unbound probability. 
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Table 2 Kinematic parameters of 19 HVS candidates. 


HVS 

d0“ 

(kpc) 

Rg'’ 

(kpc) 

Vg" 

(1cm s“^) 

Vesc-Xue'^ 
(km s“^) 

V esc —Pacy znski*^ 
(km s“^) 

Vesc—KopOSv'^ 

(km s“^) 

Veac— Kenyon‘S 
(km s“^) 

Vesc— Gnedin'* 
(km s“*^) 

1 

1.5±0.2 

8.0±0.2(0.1) 

644±96 

491 

540 

573 

592 

607 

2 

4.6±0.9 

7.5±0.9(0.3) 

626±98 

500 

540 

572 

603 

608 

3 

4.6±1.5 

9.0±1.5(0.5) 

572±195 

493 

523 

555 

595 

596 

4 

2.9±0.9 

10.1±0.9(0.3) 

563±178 

474 

519 

551 

574 

590 

5 

3.7±1.4 

8.7±1.4(0.4) 

540±218 

492 

527 

559 

593 

599 

6 

5.9±0.8 

8.7±0.8(0.5) 

524±101 

505 

524 

555 

608 

598 

7 

1.9±0.4 

7.9±0.4(0.2) 

503±83 

492 

541 

575 

594 

608 

8 

2.5±0.6 

9.4±0.6(0.3) 

501±119 

481 

524 

557 

581 

595 

9 

2.5±0.5 

8.7±0.5(0.3) 

490±98 

488 

530 

562 

589 

600 

10 

2.3±0.7 

8.5±0.7(0.3) 

489±92 

489 

532 

565 

590 

602 

11 

3.0±0.9 

9.8±0.9(0.2) 

671±115 

476 

523 

556 

576 

593 

12 

1.4±0.4 

7.9±0.4(0.1) 

621±138 

492 

541 

575 

593 

608 

13 

5.5±2.3 

13.5±2.3(0.7) 

613±157 

451 

496 

527 

548 

569 

14 

1.8±0.5 

9.0±0.5(0.2) 

603±91 

482 

530 

564 

582 

599 

15 

2.0±0.4 

8.9±0.4(0.2) 

591±60 

484 

529 

562 

585 

599 

16 

5.1±1.7 

12.7±1.7(0.4) 

583±104 

457 

499 

530 

554 

573 

17 

3.1±0.6 

8.5±0.6(0.3) 

561±93 

491 

531 

563 

592 

601 

18 

6.1±2.6 

8.4±2.6(0.8) 

527±140 

506 

526 

557 

609 

600 

19 

0.6±0.1 

8.3±0.1(0.04) 

508±42 

488 

538 

572 

589 

605 


“ Heliocentric distances obtained by the distance modulus. 

^ Galactic distances. 

Galactic total velocities. 

Escape velocities obtained by the Xue08, Pacyznski90, KoposovlO. KenyonOS and Gnedin05 Galactic 
potential model respectively. 


Beside, the unbound probability from such a Monte-Carlo simulation can only represents the probability of 
unbound to the Galaxy when kinematic parameters are in given error ranges. 


5 METALLICITY DISTRIBUTION AND POSSIBLE ORIGINS 


The metallicity distribution investigation of Galactic populations, including the Galactic bulge, the Galactic 
disk, the Galactic halo, and the globular cluster, indicate each population have a significantly different 
metal abundance distribution from others. Meanwhile, metal abundance ([Fe/H]) for a star reflects that of 
the place where it was born. So, metallicity distribution can be used as a tool to explore the origin of our 
HVS candidates. 


Sadler et al 


(Il996h measured [Fe/H] for 322 K giants of the Galactic bulge, and present the metal¬ 


licity distribution function (MDF) for the Galactic bulge. They found that the mean abundance of their 
K giant sampl es is < [Fe/H] >= —0 .11 ± 0.03, and over a half of them are in the range of —0.4 < 
[Fe/H] < 0.3. Schlesinger et al. ( 2012 ) derived the MDF of the Galactic disk using 24,270 G and 16,847 
K dwarfs from the SDSS SEGUE, different from previous investigations, this work considered obser¬ 
vational biases for the fir st time, and the ir G and K dwarf samples are the most complete samples in 


both number and volume. 


An et al 


(120131) estimated metal abundance for individual star in SDSS Stripe 
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Table 3 Probabilities that these candidates are unbound in five Galactic potential models. 


HVS 

P unbound"^Ue08 

Punbound-Pacyznski90“ 

Punbound"KopOSOV 10 

Punbound "KenyonOS 

Punbound"Gbedin05 

1 

0.93 

0.84 

0.75 

0.86 

0.64 

2 

0.80 

0.72 

0.64 

0.70 

0.54 

11 

0.87 

0.80 

0.75 

0.82 

0.67 

12 

0.80 

0.70 

0.62 

0.71 

0.53 

13 

0.81 

0.73 

0.67 

0.75 

0.58 

14 

0.80 

0.69 

0.60 

0.71 

0.51 

16 

0.77 

0.68 

0.61 

0.70 

0.51 

15 

0.89 

0.75 

0.62 

0.77 

_ b 

3 

0.61 

0.57 

0.52 

— 

— 

4 

0.64 

0.57 

0.52 

— 

— 

5 

0.57 

0.52 

— 

— 

— 

6 

0.58 

0.52 

— 

— 

— 

17 

0.69 

0.58 

— 

— 

— 

18 

0.53 

0.50 

— 

— 

— 

7 

0.54 

— 

— 

— 

— 

8 

0.54 

— 

— 

— 

— 

9 

0.51 

— 

— 

— 

— 

10 

0.50 

— 

— 

— 

— 

19 

0.63 

— 

— 

— 

— 


Notes. 

“ Unbound probability obtained by the Xue08, Pacyznski90, KoposovlO, KenyonOS and GnedinOS potential 
model respectively. 

^ ‘—’ means the candidate is bound in certain potential model, and we do not calculate unbound probability 
in this case. 


82, and presented an unbiased MDF of the Galactic halo. iHarrisI (Il997h compiled a catalogue which 


contains basic parameters of distances, velocities, metallicities, luminosities, colors, and dynamical pa¬ 
rameters for 147 globular clusters in the Milky Way, and we obtained the catalog from the website 
http://vizier.china-vo.org/viz-binWizieR?-source=VII/202fe-ref=VIZ55014f467633. 


We compare t he metallicities of our ca ndidates with the M DFs of the Gala ctic bulge (ISadler et al 


the Galactic disk (ISchlesinger et al 


20121), the Galactic halo (lAn et al 


1993, 


2013h and known globular clusters 


(lHarrislll997h . the MDFs for each Galactic population and our candidates are shown in Figure 0 From this 
figure, we can see that the metallicity distribution of our HVS candidates is well consistent with the G and K 
dwarf samples in the disk, and is also roughly consistent with the low-metallicity end of the Galactic bulge. 
Clearly, the MDF of our candidates is completely inconsistent with that of the Galactic halo and globular 
clusters. 

In addition, we estimat e Mg abundan ces [Mg/Fe] using a profile matching method in the region of Mg 
I b lines around A5170A (iLi et al.ll2014t) . and they are listed in the last column of Table [T] The external 
uncertainty of this profile match method may not be as large as 0.2 dex, and the upper limit of internal 
uncertainty is 0.3 dex estimated by the Monte-Carlo simulation. From Table [T] we can see that the Mg 
abundances of four candidates, i.e. HVS4, F1VS5, FIVS8 and HVS18, are larger than 0.5. We visually 
inspect the spectra of the four candidates, and find that the noise seriously affect spectra quality of the Mg I 
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Vgal - Vesc (km/s) 


1000 


Fig. 2 The difference distribution of V^ai and \4sc, where V"gai and \4sc are total Galactic veloci¬ 
ties and escape velocities of each sample realization in a million random samples of Monte-Carlo 
simulation. 


b region which is extremely important for the [Mg/Fe] estimation. So, we do not regard the four candidates 
as Mg-enhanced stars. For other candidates, they have high quality spectra in the Mg I b region, and their 
[Mg/Fe] fall within the range of error. Similarly, the [Mg/Fe] values of our candidates are roughly consistent 
with the Galactic bulge and disk. 

Therefore, our candidates are likely originated from the Galactic bulge and disk, and the Galactic halo 
and globular cluster may not be their possible origin place. However, the determination of exact birth place 
need much more reliable parameters, will be released in the future data release of LAMOST, to calculate 
the intersection regions of our candidates’ trajectories and the Galactic disk. 

6 DISCUSSION AND CONCLUSIONS 

In this paper, we present 19 F, G or K type hyper-velocity star candidates from over one million stars of 
the first data release of the LAMOST regular survey. We initially select over half a million F, G and K 
dwarfs with Teff and logp criteria, and then further pick out over 190,000 final F, G and K dwarf samples 
with a series of photometric criteria. Then, we obtain 6 D phase space coordinates and escape velocities 
for each dwarf, and select 17 hyper-velocity star candidates with ‘clean’ proper-motions, and 15 candidates 
with ‘reliable’ proper-motions. We finally individually inspect spectra of the 32 HVS candidates, and find 
19 of them have high quality spectra. Through checking with previous four low mass HVS catalogues in 
literatures, we conclude that they are all new findings. 
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Fig. 3 The comparison of metallicity distribution of 19 HVS candidates (black solid line) with 
the Galactic bulge (magenta), disk (green), halo (cyan) and globular clusters (blue). 


Although we use a strict criteria to ensure reliability of the kinematic parameters of our candidates, we 
still can not conhrm that we identify 19 HVS. Therefore, we calculate the unbound probability for each 
candidate using the Monte-Carlo simulation, assuming a non-Gaussian proper-motion error distribution 
and Gaussian heliocentric distance and radial velocity error distribution. Such a probability shows each of 
our candidates can escape from the Milky Way in what extent. We hnd all the candidates have unbound 
probabilities over 50%, one of them can even exceed escape velocity with over 90% probability, and the 
unbound probability varies in different potential model for each candidate. To investigate the origin of our 
candidates, we compare metallicities of our candidates with MDFs of the Galactic bulge, disk, halo and 
globular cluster, and conclude that the Galactic bulge or disk are likely the birth place of our candidates, 
and the Galactic halo and globular cluster seems not to be their possible origin places. 

When we select HVS candidates, there are a large amount of stars with dist22 < 7, which implies 
that they suffer from photometry blending from a near neighborhood. There also exist a large fraction of 
stars with nFit<5, which means that they have few position detections. The Hubble Space Telescope (HST) 
Fine Guidance Sensor (FGS) and future GAIA can provide more accurate proper-motion measurements and 
conhrm HVS from them, and they can also help us to verify proper motions for our 19 HVS candidates. 
In addition, high resolution spectroscopic observations are extremely essential to obtain more accurate 
measurements of stellar atmospheric parameters, which may help us to determine more accurate origin 
places through calculating trajectories and detailed metallicity distribution analysis, and can also decide 
whether our HVS candidates are in binaries. 

If there exist candidates in binaries, their heliocentric distances and Galactic total velocities will be 
systematically under-estimated, and their escape velocities will be correspondingly over-estimated, thus 
these binaries should more likely be able to escape from our Galaxy. Besides, we estimated the effect 
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of binary orbital velocities on the observed heliocentric radial velocities and the Galactic total velocities 
assuming three types of companion (e.g., a solar mass main sequence companion, a neutron star companion, 
or a black hole companion) in our previous work (Li, et al. 2012), and we can see that average effect of 
binary do not exceed 100 km s“^, which have little effect on our results. 
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